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We study low-energy magnetic excitations of doped spin-ladders, based 
on an effective Hamiltonian describing interactions of mobile spin and back- 



'^ ' ground spins. The helicity moduli against fluctuations are calculated in terms 

fi ' of the doping concentration, the leg number and interaction strengths in the 

^- ■ ladder. The doping ranges for spiral phase are computed and its relation to 

> 

X/-^ • spin gap structure is discussed. The spiral order can enhance existing antifer- 

^^ ' romagnetic order at certain doping range or reduce it at very low density of 

Qv ' holes. For odd-legged ladders, formation of spin gap at certain narrow doping 

-(— > 



range becomes possible via coupling of spirals and background spins. 



PACS numbers: 74.20.Mn, 75.10.Jm, 71.27.-Fa, 75.50.Ee 



Typeset using REVTgX 



There have been numerous studies on long-range magnetic order in two-dimensional 
quantum systems, some of which employed mapping to (2+l)D non-linear sigma 
model(NLsM) ||I|] for the long-wavelength, low-energy states. Adding the effects of doped 
holes to 2D antiferromagnets, several works were done of the in-plane spiral modes by mean- 
field calculation [Q ||^ |Q [Q . A couple of studies included the out-of-plane spirals generated 
in 2D systems, starting with t-J models |^ which incorporated the effects of mobile spin 
in magnetic background 0. By spiral modes, we mean vortex rotations formed by back- 
ground spin phases or long-wavelength excitations arising from fluctuations of staggered 
magnetization. In lower dimensions, some studies have been done for low-energy excitations 
in undoped spin ladders p [|1^. The ground state of even- legged ladders is known to 
be spin liquid(singlets), given strong to intermediate coupling - the exchange ratio being 
J± > -^11 • The odd-legged ladders have even-parity channels and odd-parity channel, due to 
the reflection symmetry. The mean-field theory(MFT) and numerical studies agree that the 
even-channel makes the spin liquid while the odd-channel makes the Luttinger liquid, before 
doping. After started doping, the hole spins tend to go into the even-channel, thus keeping 
the gapped mode just like even-ladders, much more likely than into the odd-channel. It 
has been mostly accepted that the "topological term" in the NLsM describing a undoped 
ladder is classified by the parity of leg numbers(n;). The term is shown to be zero for un- 



doped even-legged ladders and 2ttS rii for odd-legged ladders [|T^]. As for doped ones, there 



have been a couple of results which showed clear departure from this type of distinction, 



by DMRG calculations [|T^ and by Lanczos calculations along with MFT ||TT|. In this pa- 
per, we focus on the low-energy spin excitations in doped ladders by analytical calculations, 
discuusing the role of the topological term for doped ladders. Doped spin ladders consist 
of background spins and empty sites with residual spin. The effective interaction between 
residual spins are antiferromagnetic, if they belong to different sublattices of any bipartite 



ladder sites [|H|. New magnetic excitations arise from coupling of collective polarization and 
local magnetization, as the effective spin of the hole moves around the ladder sites. In the 
continuum limit, the distorted spin background induces non-zero "magnetic" fiux due to 



phase fluctuation. To see the spin response, the hehcity modulus is calculated as an explicit 
function of doping rate, leg number, temperature and interaction strengths in spin ladders. 
By minimizing the total energy of the system, we find the range of hole concentration which 
favor the spiral mode, from the helicity modulus. Our results of the doping range for gapped 
mode in odd- ladders agree well with those existing numerical studies [^ |]TB|], within rea- 



sonable range of interaction coupling in our model. Also, we determine how spiral modes 
enhances/suppresses AF order and the range for these effects. 

We begin with an effective Hamiltonian, based on the one made by Shraiman and Siggia 
for 2D magnets 0. At half-filling, the low energy state has the commensurate Neel order. 
For the long-wavelength staggered spin state, the Neel vector is twisted, by aj^ around z-axis 
at the end of the ladder length L with respect to the first n and by ^7;;-r[vr ^^ ^^^ ^^^ °^ 
the ladder with n^-leggs (with the lattice spacing oq). When the mobile spins are introduced 
upon doping, the Neel vector is allowed to rotate around y-axis as well, to derive all modes 
generated by collective polarization and phase fluctuation, n = (sin 6 cos (p, sin 6 sin 0, cos 6) . 
This is equivalent to twisting the spin quantization axis from site to site via appropriate 
transformation as done in some existing works [ffl P| H]- The effective Hamiltonian is 
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where the magnetization current ja{(l) = n x Van and the dipole polarization Pa{q) = 
X^fcSin A;a^;j,_2 f'^^,^,. The second and the third term is the rotational kinetic energy of 
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the unimodular n. The subscript a stands for in-plane(IP) mode and out-of-plane(OP) 
spirals. The wavevector of the mobile spin is related to the spiral wavenumber Qa as 
< ^k^2 \t';j, 2 >^ QaSinka- Identifying the components of incommensurability wavevec- 
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tor Q = {Qxi Qy) with the angular terms in the effective interaction term in H^ff, the spiral 
wavevector is Q —>■ — sin^ cos^(cos0 + sin</))] = cos{Qx/2) — cos{Qy/2). The spiral wavevec- 
tor has a symmetry similar to a d-wave. The effective propagation in real space is along the 
rung direction of the ladder, having arisen from triplet formation along the ladder direction. 



This mode is equivalent to the transverse mode |^ or to the stripe phase |^ as termed for 
2D systems . The doping range for the spiral order to occur is computed from the helicity 
modulus against IP twisting. The spiral wavenumber Q is proportional to the doping rate 
0. For the purpose of calculating the helicity modulus, we use Q ~ 5y- for IP mode 
and Q ~ (5 for OP spirals. For IP mode, the effective Hamiltoinian density is given 
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The free energy per unit area of the ladder plane is 
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where A = —2g6j^, B = —2g/6j=- and (3 = l/fc^T. The subscript ± refers to those along 
the rung direction while || is for those in the chain direction. The spin wave velocity for 
doped ladders is derived from mapping the unit ladder cell (repeated regularly as defined by 
the density of empty sites) to NLsM. Taking only the slow component of the Neel vector and 

for the coherence length ^ the ladder width(aon;), the coupling now includes the doping 

1 

effect g = [{1 — 5)r7,;5']^^ [^ ^ 2j^] ■ "^^^ ^P^^ wave velocity is collected from the prefactor, 
Vs = 2SJ\\ ("tix ) ( 1 + 9^ ) • As for 2D systems without any hedgehogs, it is known that the 
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topological term — ^ / [damYcPr is zero, regardless the spin number or the configurations 
of lattice sites |15| ||16[ [1^ . We discuss the different meaning of this term in doped ladders 
later. The coupling in the effective interaction(V^/) is found from relating the hole density 
to the collective polarization < Pg, > and the incommensurability wavenumber Q of the 
spiral mode. From < Vg, >^ 2gf Qln[2(n/ — 1)] < P >~ gf Q{ni — l)L/167r^ao, we estimate 
the coupling gl ^ J_lQ/ < P > ~ 32ao7r^ ln[2(?7,; — 1)](5J_L/L(n/ — 1), where L is the length 



of unit ladder cell, which decreases as hole density increases. The other coupling g can be 
estimated from renormalized coupling((7o-) of the NLsM; 1/g = l/g„+TT[l/{k'^ + V{g)}]. For 
instance, given the most realistic interaction strengths in cuprate ladders j^ = 5,g = 1 with 
S = 0.1, while 5 ^ 1 for much greater g. The helicity modulus(HM) against the in-plane 
fluctuations is computed 
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This is an explicit function of the doping concentration 5, leg number n^, temperature and 
the interaction strengths in the ladder t^, J±_ and J\\. Fig. 2 is a surface plot of the IP 
HM versus doping rate and temperature. The common features of HM for all ladders are 
(i) that the HM diverges at certain doping rate 5 = l/{2gj^), (ii) that the sudden jump is 
followed by sharp decrease in HM to negative values, (iii) that the fiat plateau is over wide 
range of doping and temperature, (iv) that another range for shrap increase of HM is at 
0.8 < S < 0.9. The spin liquid ground state acquires excited triplets upon doping. When the 
hole density reaches a certain level, intrachain pair formation becomes substantial as well, 
which makes the whole system very stiff. The diverging modulus followed by drastic decrease 
as well as flattening over wide range, indicate presence of spin gap, since the HM(stiffness) 
goes as the inverse of spin susceptibility. In the region of negative HM, the system releases 
spiral modes to lower the energy. Though HM itself is a local response, the spiral phase 
(bulkwise) can be formed in the regime of negative HM. As doping increases further, there 
is another sharp increase in HM at very low temperatures, due to coupling of spirals and 
intrachain spins. The spiral modes can put a pair of spins into a level in the gap. Thus, 
the short-range AF order of spin liquids can be lost at certain doping regime. In even- 
ladders with existing spin gaps, spiral modes can reduce the spin gap sizes. At higher g 
which is also physically feasible, the spiral phase is absent at low doping, though the other 
window at much higher doping still remains. In Fig. 2 for 3-leg ladder, the possibility of 



generating spin gap in odd-ladders exists around S = 0.1 with g=l. For more attractive 
g > I which is very hkely in physical ladders, the spiral mode and gap formation occur 
in smaller doping regime oi 6 < 0.1, at temperatures T < 5K. The region where the HM 
becomes negative is commonly --^ < 5 < 0.5 at T < 5K, producing spiral phases for 
all ladders. The characteristics are repeated for ladders with any number of leggs, though 
precise values are a bit different, due to the leg-number dependence in g. As leg-number 
increases, the doping regime for IP spiral mode is reduced. Since the even-parity channels 
in odd-ladders are dominant over the odd-parity channels, the added (effective) spins tend 
to be in the even-channels. Given the attractive potentials provided by V^'s in our model, 
binding of the spins in the odd-channel and the spirals(solitons) yield an energy gap. This 
gives rise to a gapped phase at 6 > Sspirai in odd-ladders, with all channels participating. 
Our result of doping rate for the gapped mode in odd-ladders {Sspirai = 0.1 for ^ = 1) is close 
to the result by Rice et. al. |jll| who found the 6c > 0.13. If the coupling ^ > 1 in physical 



ladders, which gives our 6spirai < 0.1, the DMRG result of 6c > 0.06 by White et. al |[T^ is 
comparable as well. In our picture, there is interplay between effective interaction strengths 
and hole concentration, giving the coupling g vary consequently. Existing interpretation 
of gap generation in odd-ladders by White et.al. uses the idea of domain walls formed by 
competition between the kinetic energy and the exchange energy. In our picture, domain 
walls can be interpreted as areas of the spiral vortices with different spiral wavenumber or 
different chiralities. Alternative understanding has been given with proximity effect to the 



gapped mode of the even-channel, which enhances the pairing within the odd-channel ||Tl 

For out-of-plane(OP) mode, the effective Hamiltonian for OP spiral is Hqp 
{da'ip)'^ ± 2g/6sinkcoskcos6. The free energy per unit area along the z is F 
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For a exchange ratio of J± = 2J\\, the result shows formation of OP spiral phase over wide 
doping range < 6 < 0.5 for ladders with any number of rii > 1. In all ladders, there are 
sharp increases in OP HM at 6 > 0.5 for T < 15K, making OP spirals robust. At very 
high hole density beyond this 6, the peaked modulus indicates forming bound states of OP 
spirals with opposite chirality(soliton-antisoliton pairs). The spin quantum number would 
be integers for this new type of condensate. The HM(stiffness) decreases as temperature 
increases at high doping regime for all ladders. From minimizing the Hamiltonian H^ff with 
respect to a mean- field Q [0], the self-consistency condition gives Q = ^ J2k sin k {n'^ + n~) ~ 
^ J2k sin k {^ — Sk)- The Fermi distribution function n^ = exp[/3(£:fc =F ^ Q sin fc — /x)] -|- 1 
gives 71^ — n^ k, 2f3g Q sin k at low T. Then the mean effective interaction terms are 
<V~> + <Vg, >~ -% Efc sin^ k {2/3g Q sin k)'^ + ^^Q'^. The < Hmf > becomes positive 
when the helicity modulus 3? becomes negative, which makes the system unstable. So, the 
system releases the torsional momenta. That is, the spiral phase with uniform wavenumber 
throughout the system is favoured when 3fJ < 0. From 3fJ/p, we find the doping range for 
spiral phase --^ < 5 < |, which agree with the graphical result. The spin coherence length 
in the ladder plane is computed when the HM is zero, 

128ao7rln[2(n,-l)]^/3J^ 

When the incommensurate spin state is stable in the regime of positive HM, it is possible 
to gather the (hole) spins along valleys of the spiral wavevector. This may result in phase 
separation of hole-rich region and no- hole region [Q |1^, adjacent to the spiral phase. As for 



spin ladders with even number of leggs, the upper limit of doping for Luttinger liquid phase 
is ~ J/2tg. The lower limit would be around J/4:tg. In our results, there is no transition 
like Kosterlitz-Thouless(KT) transition in 2D ||20|, since there is no strong size dependence 
and no plateau in the stiffness around any particular temperatures. Though there are size 
dependent terms in the expression of the HM, these are rather weak. Calculated spin 
coherence suggests that the gap function has a different form from this, due to the doping 
effect. As far as the helicity modulus, ladders with even/odd-number of leggs do not have 



differing signs, since the factor ln[2(?2i — 1)] does not distinguish the parity in total leg number 
111. This is the difference from the result on ID spin chains which distinguish the parity, but 
analogous to 2D systems which do not differentiate the spin quantum number being integer 
or half-odd integer. For ladders with finite number of leggs, the boundary condition in the 
spin wavefunction is given \E';(?2 + 1) = ± exp[— i0]\l'i(l), where / is the leg index and n is the 



site number in the block of the system [^. The meaningful parity in ladders is that of (leg 
number plus site number), instead of the leg munber alone. The "+" or "-" sign goes as the 
chirality associated with the winding orientations of spirals, having arisen from the presence 
of two sublattices(of up or down-spin). So, the ground state in odd — ladders is degenerated 
for spiral modes with different chiralities at different ladder regions, thus satisfying Lieb- 
Mattis-Shulz theorem. This means that gapped mode is possible in odd-ladders. This is the 
same result as found by previous MFT study on odd-ladders [|14| . 



So far, we restricted the spiral winding number (w) to 1, to ensure the spinor wavefunction 
to be periodic (instead of anti-periodic). Corresponding staggered spin flux(chirality) in 
continuum limit is 2itw. In terms of the Neel vector at bipartite ladder sites, the winding 
number is = [ui ■ n2 >^ n^ + fli ■ n^ x n4 — fii ■ n2 >< n^ — n2 ■ n^ x n^]. If the systems were 
undoped ladders described by NLsM, the w would be equal to the topological term. This 
gives the solid angle subtended by the Neel vectors or the Berry phase acquired by the 
Neel vector while traversing around the path in spin space ||22| ||2^. Allowing the spiral 
modes, the winding number need not be integers. This implies that the topological term in 
doped ladders does not simply determine their classes by the leg numbers. The induced field 
acting as twisted boundary would affect the magnetic ordering in ladders. To see whether 
spiral modes compete with antiferromagnetism in ladders or not, we rewrite the Neel vector 
including the winding number, n = {sm6cos{w(j)),sm6 sm{w(f)),cos6). Then, the effective 
interaction term is Vr ~ —n'^wgS, coming from the gradient of n. This can provide attractive 
potential |2^ between the polarised spin and the spin current if the w is positive. To maintain 
the attractive potential between the spin texture and the spiral (soliton), w is chosen to be 
of the chirality "-I-". This would result in destroying the AF order of background spins. 
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However, there is another interaction term Vg, in our effective Hamiltonian, which does not 
depend on the winding number. This contribution can be greater in length scales of a few 
lattice constant, though it decreases as the doping decreases and/or the leg number increases. 
In cases of very low hole density (a few percent) and many number of leggs, it is possible 
that the V~ wins to reduce the AF order. A known result by Arrigoni et. al. on spiral order 
in 2D Hubbard model shows that the transverse mode enhances AF order. The other spiral 
mode can enhance ferromagnetic order of the background. Another work by Mori et. al. 
on 2D t — J model also show two similar spiral phases. As for doped ladders which started 
with commensurate Neel order at half-filling, our result show that the whether spiral modes 
reduce or enhance the AF order depends directly on the hole density. Only at very low 
doping regime (a few percent), spiral order suppresses antiferromagnetism. It is due to that 
the effective interactions Vg/ and V~ compete in some doping regime. 

In summary, due to the hole motion distorting the spin texture in doped ladders, phase 
fluctuations of the staggered magnetization yields new type of excitations. In continuum 
limit, long-range spiral modes arise and give the helicity modulus vary depending on the 
doping rate, temperature, leg number and interaction strengths. The spiral modes can put 
magnons into sublevel(s) in the gap. The spin gap is reduced in even-ladders while it can 
be generated in odd-ladders, depending on hole density. Since the even-parity channels in 
odd-ladders are dominant over the odd-parity channels, the added spins tend to be in the 
even-channels. Thus, in the begining of doping, the even-channels give the spin gap as 
spin liquids of undoped ladders. Given the attractive potentials provided in our model, the 
spin current in the odd-channels tend to make bound states with the spirals in the even- 
channel. The enhanced pairing brings out the gapped Luther-Emery phase at 5 > Sspirai 
in odd-ladders. Our results of the doping range for gapped mode in odd-ladders agree well 
with those existing numerical studies. Antiferromagnetism in ladders can be enhanced over 
substantial doping range via having spiral phase. At very low doping of few percentage, 
spiral modes promote ferromagnetism. In this paper, we showed analytically, generation of 
spin gaps in odd-ladders and the explicit relation between magnetic orders and topological 



excitation. In our forthcoming study, we are to examine whether the spiral phase enhances 
superconducting transition or suppresses via tunnehng of sohtons. 
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FIGURES 




FIG. 1. Spiral wavevector Qy vs. Q^ for < < vr. Q reaches the maximum at Qx = 0, which 
indicates the gap formation along the y-direction. 
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FIG. 2. IP hehcity modulus(HM) for 3-leg ladder with g = l,t^ = 0.2 eV and J_l = 0.04 eV. 
The temperature is given in K. When HM becomes negative, spiral phase is formed at 0.1 < 5 < 0.2. 
The sharp increase at higher 5 > 0.8 with subsequent negative region indicates another window of 
spiral mode. As leg-number increases, the doping regime for IP spirals is shifted to smaller values. 
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FIG. 3. OP helicity modulus for 3-leg ladder. OP spiral mode is formed at doping range 
Q < 5 < 0.5. Overall features are common for ladders with more number of leggs, given the weaker 
dependence of the HM on leg-number than the dependence on doping or temperature. 
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